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gauge transformations and gauge invariance 



m 
o 
o 

(N 



c3 



o 

O 



o 
cr 

wo. 



X 



Marco Bruni| Leonardo Gualtieri§ and Carlos F. SopuertaJ 

X Institute of Cosmology and Gravitation, Mercantile House, University of 
Portsmouth, Portsmouth POl 2EG, Britain 

§ Dipartimento di Fisica "G. Marconi", Universita di Roma "La Sapienza" and 
Sezione INFN ROMA 1, piazzale Aldo Moro 2, 1-00185 Roma, Italy 



Abstract. An implicit fundamental assumption in relativistic perturbation 
theory is that there exists a parametric family of spacetimes that can be 
Taylor expanded around a background. The choice of the latter is crucial to 
t"""""" ' obtain a manageable theory, so that it is sometime convenient to construct 

£Nj ' a perturbative formalism based on two (or more) parameters. The study of 

perturbations of rotating stars is a good example: in this case one can treat 
, the stationary axisymmetric star using a slow rotation approximation (expansion 

■ in the angular velocity Q), so that the background is spherical. Generic 

|^r^ ■ perturbations of the rotating star (say parametrized by A) are then built on top 

f*-) of the axisymmetric perturbations in Q. Clearly, any interesting physics requires 

non-linear perturbations, as at least terms AH need to be considered. In this 
paper we analyse the gauge dependence of non-linear perturbations depending 
on two parameters, derive explicit higher order gauge transformation rules, and 
define gauge invariance. The formalism is completely general and can be used in 
different applications of general relativity or any other spacetime theory. 



PACS numbers: 04.25. Nx, 95.30.Sf, 02.40.-k 
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1. Introduction 

An implicit fundamental assumption in relativistic perturbation theory is that there 
exists a parametric family of spacetimes such that the perturbative formalism is built 
as a Taylor expansion of this family around a background. The perturbations are then 
defined as the derivative terms of this series, evaluated on this background P . In most 
cases of interest one deals with an expansion in a single parameter, which can either 
be a formal one, as in cosmology [21 El El or in the study of quasi-normal modes of 
stars and black holes [HIE], or can have a specific physical meaning, as in the study of 
binary black hole mergers via the close limit approximation [7|E], or in the study of 
quasi-normal mode excitation by a physical source (see E] and references therein) . 
Typically the perturbative expansion stops at the first order, but recent interesting 
developments deal with second order perturbations El El E3 \U\ ■ 

In some physical applications it may be instead convenient to construct a 
perturbative formalism based on two (or more) parameters, because the choice of 
background is crucial in having a manageable theory. The study of perturbations 
of stationary axisymmetric rotating stars (see |12l 1131 H^] and references therein) is 
a good example. In this case, an analytic stationary axisymmetric solution is not 
known, at least for reasonably interesting equations of state. A common procedure is 
to treat axisymmetric stars using the so-called slow rotation approximation, so that 
the background is a star with spherical symmetry |151 116j . In this approach the first 
order in fl discloses frame dragging effects, with the star actually remaining spherical; 
ft 2 terms carry the effects of rotation on the fluid. This is intuitive from a Newtonian 
point of view, as rotational kinetic energy goes like fl 2 . This approximation is valid 
for angular velocities fl much smaller than the mass shedding limit flu = y M/Rg tar , 
with typical values for neutron stars fix ~ 10 3 iJz. Therefore the slow rotation 
approximation, despite the name, can still be valid for large angular velocities. In 
practice, the perturbative approach up to fl 2 is accurate for most astrophysical 
situations, with the exception of newly born neutron stars (see |17| and references 
therein) . 

Given that the differential operators appearing in the perturbative treatment of 
a problem are those defined on the background, the theory is considerably simplified 
when the latter is spherical. Generic time dependent perturbations of the rotating star 
(parametrized by a dummy parameter A and describing oscillations) are then built on 
top of the stationary axisymmetric perturbations in fl. Clearly, in this approach 
any interesting physics requires non-linear perturbations, as at least terms of order 
Xfl need to be considered. A similar approach could be used to study perturbations 
of the slowly rotating collapse, even if in specific cases ^3 ^3 E0] the perturbative 
expansion depends by one parameter only. 

Classical studies in the literature have not analysed in full the gauge dependence 
and gauge invariance of the non-linear perturbation theory. For example, in [IBj 
the second order perturbations are treated in a gauge invariant fashion on top of 
the first order perturbation in a given specific gauge. The perturbation variables 
used are therefore non gauge invariant under a complete second order gauge 
transformation |21 12 1|. but only invariant under "first order transformations acting 
at second order" JHj- While this may be perfectly satisfactory from the point of 
view of obtaining physical results, one may wish to convert results in a given gauge 
to a different one El> to compare results obtained in two different gauges, or to 
construct a fully gauge invariant formalism. To this end one needs to know the gauge 
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transformation rules and the rules for gauge invariance, either up to order n |2] or at 
order n only, as in The situation is going to be more complicated in the case of 
two parameters, as we shall see.J 

In this paper we keep in mind the above practical examples, but we do not 
make any specific assumption on the background spacetime and the two-parameter 
family it belongs to. As in El EEj, we do not even need to assume that the 
background is a solution of Einstein's field equations: the formalism is completely 
general and can be applied to any spacetime theory. We analyse the gauge dependence 
of perturbations in the case when they depend on two parameters, A and Q, derive 
explicit gauge transformation rules up to fourth order, i.e. including any term X k il k 
with k + k' < 4, and define gauge invariance. This choice of keeping fixed the total 
perturbative order is due to the generality of our approach. In practical applications 
one would be guided by the physical characteristics of the problem in deciding where 
to truncate the perturbative expansion. For example, in the case of a rotating star 
one could consider first order oscillations, parametrized by A, on top of a stationary 
axisymmetric background described up to fl 2 , neglecting therefore X 2 fl terms. Or 
instead, one could decide that A 2 S1 terms are more interesting than the Af2 2 ones in 
certain cases. 

From a practical point of view, our aim is to derive the effects of gauge 
transformations on tensor fields T up to order k + k' = 4. It is indeed reasonable 
to assume that in a practical example like that of rotating stars, at most one will want 
to consider second order oscillations ~ A 2 on top of a slowly rotating background 
described up to (3(f2 2 ), in order to take into account large oscillations and fluid 
deformations due to rotation. 

We will show that the coordinate form of a two-parameter gauge transformation 
can be represented by: 

* p = x>i + ^£(1,0) + ^£(o,i) 

+ ~2~ (£(2,0) + £(l,0)f(l,0),L/) + ~Y (^(0,2) +£(04)£(04)^) 

+ ™ + eo^i,o)£fo,i)„ + ^(0,1)^,0), J + 3 (A, n) , (1) 

where the full expression is given in Eq. (|75l) . Here x y £K ^, and 

£(02) are independent vector fields and (eo,ei) are any two real numbers satisfying 
go + ei = 1. Coupling terms like the Xfl in (JIJ are the expected new features of the 
two-parameter case, cf. [21 124| . Our main results are the explicit transformation rules 
for the perturbations of a tensor field T and the conditions for the gauge invariance 
of these perturbations. 

The paper is organized as follows: in Section |2 we develop the necessary 
mathematical tools, deriving Taylor expansion formulae for two-parameter groups 
of diffcomorphisms and for general two-parameter families of diffcomorphisms. In 
Section|21we set up an appropriate geometrical description of the gauge dependence of 
perturbations in the specific case of two-parameter families of spacetimes. In Section^] 
we apply the tools developed in Section El to the framework introduced in Section |21 
in order to define gauge invariance and formulas for gauge transformations, up to 

X The concept of perturbation theory with more than one parameter has already been introduced, 
for first-order perturbations, in 1221 . where the standard definition of spacetime perturbations 1231 is 
extended by using a (4+n)-dimensional fiat space in which space-times are embedded. The main aim 
of these works is to re-examine the gauge invariance of the metric. 
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fourth order in the two-parameter perturbative expansion. Section [S] is devoted to 
the conclusions. We follow the notation used previously in [21 1241 125) for the case of 
one parameter perturbations. 



2. Taylor expansion of tensor fields 

In order to consider the issues of gauge transformations and gauge invariance in 
two-parameter perturbation theory we need first to introduce some mathematical 
tools concerning the two-parameter Taylor expansion of tensor fields. Since Taylor 
expansions are aimed to provide the value of a quantity at some point in terms of its 
value, and the value of its derivatives, at another point, a Taylor expansion of tensorial 
quantities can only be defined through a mapping between tensors at different points 
of the manifold under consideration. In this section we consider the cases where such 
a mapping is given by a two-parameter family of diffeomorphisms of M. , starting from 
the simplest case in which such a family constitutes a group. 



of 



(2) 



2.1. Two-parameter groups of diffeomorphisms 

Given a differentiable manifold M, a two-parameter group of diffeomorphisms 
M can be represented as follows 

<t> : M x H 2 — > M x H 2 
(p,A,n) (0x,n(p),A,n). 

For the purpose of introducing two-parameter perturbation theory, where perturbing 
first with respect to the parameter A and afterwards with respect to the parameter f2 
should be equivalent to the converse operation, we will assume that tpx.n is such that 
it satisfies the following property§ 

0Ai,Oi ° 0a 2 ,o 2 = 0Ai+a 2 ,Oi+o 2 > V A, G H . (3) 
By one hand, this property implies that the two-parameter group is Abelian. On the 
other hand, it allow us to make the following useful decomposition of 0a, n into two 
one-parameter groups of diffeomorphisms (flows) that remain implicitly defined by 
the equalities 

4>x.n = 0A,O ° 00,0 = 4>0. Q ° 0A,O ■ (4) 

The action of the flows (j)\,o an( i 00, fi is generated by two vector fields, r\ and £ 
respectively, acting on the tangent space of M x K 2 . The Lie derivatives of a generic 
tensor T with respect to r\ and £ are 



C V T 



= lim — (i 

A^O A 

= lim i( 

f2->0 £1 



y A,0 



T-T) = 



J o,n 



T-T) = 



d_ 

d\ 

d 

dJ} 



T 



—<l>on T 





' d 


. A=0 





0=0 



d_ 



(5) 
(6) 



where the superscript * denotes the pull-back map associated with the corresponding 
diffeomorphism 0. Because the group is Abelian the vector fields r] and £ must 
commute 



(7) 



§ We adopt the convention that the first label in <j> corresponds to the flow generated by r\ and 
parametrized by A, while the second label corresponds to the flow generated by £ and parametrized 
by U. Therefore <f> x ,n ^ (f>n,\- 
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The Taylor expansion of the pull-backs cj>* x Q T, 4>q n T is given by (see [2]) 



00 X k 
■o T = Elf 



k=0 



k=0 



n k 

~kl 



dX k 



J X,0 



T 



° x k 



IJ o,n 



T 



A=0 



V — C k T 



fc=0 

00 fe 



E 



(8) 



(9) 



From this, using we can derive the Taylor expansion of the two-parameter group 
of pull-backs <^> A qT: 

X k Q k 



E 



fc,fc'=0 



dxkmW vx ' Q 



T 



00 



k,k'=0 



klk'l " c 



(10) 



12. Two-parameter families of diffeomorphisms 

Let us now consider the general case of a two-parameter family of diffeomorphisms $: 



$ : X x nr 

(p,A,n) 



M x R 
(*A,n(p),A,fi). 



(11) 



In this generic case the diffeomorphisms $a,o do not form in general a group. In 
particular, 

which means that we cannot decompose $a,o as in the case of a two-parameter group 
of diffeomorphisms. The Taylor expansion of the pull-back of <&a n is formally given 

by 

Qk+k' 



fc,fc'=0 



A fc fi fc 
fcffc 7 ! 



(12) 



Since the diffeomorphisms $a,o do not form a group we cannot write this expansion 
directly in terms of Lie derivatives as in the previous case, Eq. l|10fl . Nevertheless, 
in order to study the characteristics and properties of the gauge transformations 
to be derived in Section 0] we would like to find an alternative way of expressing 
the expansion (|12fl in terms of suitable Lie derivatives, in a similar way as it was 
done in [21 US] in the one-parameter case. To this end, new objects called knight 
diffeomorphisms were introduced in |2] I25| . Broadly speaking, they are a composition 
of one-parameter groups of diffeomorphisms that can reproduce the action of a given 
family of diffeomorphisms, in such a way that the more groups we compose the better 
the approximation is. When the number of composed groups tends to infinity the 
reproduction of the family is exact (see [3 1241 05] for a detailed account on knight 
diffeomorphisms). The main aim for using these objects was to show [5]|21] that it is 
possible to expand a one-parameter family of diffcomorphism, at orden n, using Lie 
derivatives with respect to a finite number n of vector fields. Knight diffeomorphisms 
constitute an elegant formulation of this kind of expansions. In order to apply these 
ideas to our case, we can think of translating the idea of knight diffeomorphims from 
the one-parameter to the two-parameter case. We have studied this question and we 
have found that in the two-parameter case there are several ways of formulating knight 
diffeomorphisms, and they lead to different formal expansions of a two-parameter 
family of diffeomorphisms. Actually, we have found that some of these formulations 
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can be inconsistent in the sense that the differential operators that come out from them 
do not satisfy the Leibnitz rule. Then, we cannot associated to any of these operators 
a vector field whose Lie derivative will describe the action of that operator. Despite 
of these facts, we find that the main goal for introducing knight diffeomorphisms can 
be achieved in the two-parameter case. That is, we can still expand a two-parameter 
family of diffeomorphisms, at order n, using Lie derivatives with respect to a finite 
number of vector fields. 

Then, here we are going to show how to approximate a given family of 
diffeomorphisms $, up to order X k il k with k + k' = 4, in terms of some differential 
operators that we will show later can be identified with Lie derivative operators with 
respect to some vector field. In other words, we are looking for an expansion of the 
following type: 



A fe f2 fe (Combination of differential operators)^ k , T . (13) 

k,k' = 
k + k' = 4 

We have performed this expansion term by term. For the sake of brevity we present 
here only the most relevant developments and results. We start by introducing the 
set of differential operators that we have used to build such an expansion by giving 
their actions on a general tensorial quantity T:\\ 

d 



£ (2,1) T 



d 2 



d 2 



dXdfl 
8 2 



Q3 



Q3 

— $t oT 



r 2 t 

' ( e o-C(i,o)£(o : i) + ei^(o,i)^(i,o)) T . 

' ^(0,1) ^ ' 
3£(1, 0)^(2, 0)^ ~ £(1,0)^ ' 



" 2£(i,o)^(i,i)2 n — £(0,1)^(2,0)^ ~ 2e 2 £(i i0 )£(o,i)i3(i i o)j 
- (ei - e2)£(o,i)£(i,o) T - ( £ o - e 2 )£( 1)0 )'C(o ! i)r, 

Q3 

axon 2 A ' a 



(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



(20) 



~~ 2£(o, i)£(i,i)T — £( 1)0 )^(o, 2)T — 2e 3 £(o : i)£(i, o)£(o,i)^ 

- (e - e 3 )£(i )0 )£( 0) i)T - (ei - e 3 )£( 0)1) £(i )0 )T , (21) 

|| The subscripts (p, q) denote the lowest order in the expansion 1131 at which these operators will 
appear for the first time. See equation I3UI below. 
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£(4,0)^ 



£ (3,1) T 



£(1.3)T 



£(0,4)?" 



3£(o, 1) Ao, 2)T — £( j)T . 



6 £(1,0)£(2,0) T ' 



L (i,o) 



T. 



7 

(22) 
(23) 



<9 4 



d\ 3 dn x ' n 

- 3£(i,o)£(2, i)T — C( 0t i)C(3$)T — 3e4£( 2: o)£(i,i)r — ^5^-{i,i)^{2,o)T 
' 3^(1,0)^(1, l)^ — 3 (£o£(i :0 )£(o,i)^ + £ i^(o, i)£(i,o)) £(2,o)T 

- (ei - £2 — £6)£(o, i)£(i, o)^ 1 ~ 3e 6 £(i :0 )£(o, i)£(i o)2" 



- 3(e 2 - e6)£(i ; o)£(o,i)£(i,o)^ 1 _ ( e o — 2e 2 + £6)£(i )£(o,i) , 
d 4 



(24) 



_ 2£(1,0)£(1,2)-' 



T + 2£( 0il )£( 2i i)T - 2£ 



2 



T 



~ e 7<£(2, 0)^(0, 2)T — £8^(0, 2)^(2, 0)^ — ^(1,0) ^(0,2) ^ — £( , 1) ^(2,0) ^ 

- 4 (e £(i : o)£(o,i)r + £i^(o,i)£(i,o)) £(1,1)^ 

+ (£3 + £2 - £l + £9)^(0, 1)^(1, 0) T + ( £ 3 + £2 - £0 - e 9)£(i,o)£(0,l) T1 

- 2(e 3 + £ 2 - £o£l - £g)£(i i0 )£(0, 1)^(1,0)^(0, 1)^ 

- 2(e 3 + £2 - £o£i + £9)^(0, 1)^(1,0)^(0, i)£(i,o)^ 

+ 2(f 3 - £ £i)£(i i0 )£(o, i)£(i, o)T + 2(e 2 - £ £i)£(o,i)£(i, )'C(o,i)? 1 > (25) 
d 4 

- 3£(0,1)£{1,2)T — ^(1, 0)^(0,3)^ — 3£i £(o,2)£(l,l) T — Sfin^fx^j-Cjo^)? 1 

_ 3£( , 1)^(1, i)2" ~ 3 (£o£(i,o)£(o,i)^ + £i^(o,i)£(i,o)) £(o,2)T 

- (£0 - £3 - £l2)£(l, 0)^(0,1)^ _ 3£i 2 £(o,l)£(l, 0)^(0,1)^ 

- 3(£ 3 - £i2)-C 2 01) £(i,o)£(0,l)7 1 - ( e l - 2£3 + £12)^(0, 1)^(1,0) 

d 4 

(S* j 1 

on 4 x n \ 

- 4£( 0il )£( 0:3 )T - 3£( 2) T - 6£( 01) £(0,2)7" - 



(26) 



'(0,1) 



where the quantities ca {A 
following conditions 



(27) 

0, . . . , 12) are real numbers which must satisfy the 



£0 + £1 = 1 , £4 + £5 



1 



£7 + £8 



£10 + £11 = 1 



(28) 



These constants reflect the fact that, as already mentioned above, there is not a unique 
way of constructing an expansion of the type l|13|) , actually the number of possibilities 
is infinite. In this paper we have restricted ourselves to expansions in which every 
single term has the form £( pi gi ) • ■ • C( Pn <qn )T with Pi + <Zi < ■ ■ • < p n + q n . Then, 
the constants €a express the freedom that we have in constructing the expansion (|13|l 
with this criteria. This freedom and the non-uniqueness of the construction above is 
not a problem: the operators £( Pl9 ) and the corresponding vectors (see below) that we 
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are going to define should be though of as a "basis" for the construction in terms of 
Lie derivative of each of the terms in the expansion (|12|) . This "basis" is not unique, 
but the result at each order is. 

On the other hand, it is straightforward to check that these operators are linear 
and satisfy the Leibnitz rule and hence they are derivatives. We can also check that 
they satisfy the rest of conditions of the theorem stated in |Appendix A| Therefore, 
for each of them there is a vector field £(p,g) such that 

£^ T --= C M T (P^GIN). (29) 
In the particular case when $ is a group of diffeomorphisms we recover the 

previous case (subsection 12 and £( p . g ) = if p + q > 1. 

Using the differential operators we have just introduced we can express the Taylor 

expansion (|12fl of <&a.si, up to fourth order in A and Q, in terms of the Lie derivatives 

associated with the vector fields l|29l) : 

n,nT = T + \C 6(m T + nC i(0A) T 

+ y { A (2 , 0) + 4 ( i, n > } T + t + 4o.i) } T 

+ Xn + eo/^ci.o^to.i) + ei£ S(01) £ S(10) } T 

A 3 r i 

+ y \ £ £(3,0) + 3£ £(1,0) £ £(2,0) + 4(1,0) } T 

~\ 2 i + £(i,o) "^£(1,1) + £(o,i) -^£(2,0) + 2 e 2^^ (li0) £^ (01) £^ (1 0) 

+ (£1 - e 2 )A(o,i> 4(1,0) + ( £ ° - ^Kfa.oAo.i) } T 

AO 2 

"I 2~ {^"£(1,2) + ^£(o,i) "^£(1,1) + 4i.o> -^£(0,2) + ^ e 3-^$(o,i) -^--€(1,0) ^£(o,i) 

+ (e - e 3 )A ( i .o)4o ,d + ( £l " e 3) £ f(o.i) £ £(i,o) } T 
SI 3 r 3 -i 

-Q \ £ £(0.3) + 3£ £(0,1) ^£(0,2) + 4(0,1) / T 

{ £ ?(4, , +4^(1,0)^(3,0) +3 4(2,0) + 6 4(1,0) £ £(2.0) + 4(1.0) } T 



G 

A 4 
+ 24 
A 3 f2 

"I g~~ { -4.3,1) + 3 ^£(1,0)^£(2,1) + -^£(0,1) ^£(3,0) + 3e 4^5 (2 , 0)-^C(l,l) 

+ 3e 5 £ S(11) £ S(2i0) + 3£ 2 (i o) £ 4(1 1; + 3 (e £f (1 , 0) £{ (0 1} + exC^ 0:1) C^ 10) ) C^ 2 0) 

+ (e X - £ 2 - e 6 )A(o,i)4(i,o) + 3e 6 £ £(i,o) A(o.i) 4(i,o) 

+3(e 2 - e 6 )£ 2 (i o) /: C(0 0) + (e - 2e 2 + e 6 )4 (1 , 0) A(o ,d } T 



A 2 ft 2 



{■^£(2,2) + ^£(1,0) ^£(1,2) "I" ^£(0,1) ^£(2,1) + ^4(1,1) 



+ ^4(2, 0) 4o,2) + e 8^C(o,2) ^4(2,0) + 4(1,0) 4(0,2) + 4(0,1) 4(2,0) 

+ 4 (e £ 5(1 _ 0) £ 5(01) + ei£ ?(0;1) £ ?(10) ) a) 

-(es + e 2 - ei + e 9 )£ 2 (o i) £ 2 (i _ o) - (e 3 + e 2 - e - e 9 )4 (1 , ,4 ( o 

+ 2(e 3 + e 2 - eoei - e 9)A ( i, )- C £(o,i)- C £(i, 0)4(0,1) 

+ 2(e 3 + e 2 - e ei + £9)4(o, i)4(i,o)4(o,i)4(i,o) 
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-2(e 3 - e ei)C i{10) Cl (o i) C i(l g) - 2(e 2 - e Q e x )C^ A) C\ {1 0) C i{0 1} | T 
Ail 3 

^~ { £ £(1,3) "I" ^f(O.l) £ f (1,2) £ £(1,0) £ £(0,3) + ^ e 10^4(O,2) ^4(1,1) 



+ 3en£ ?(11) £ ?(0 2) + 3£| o i) £ ?(11) + 3 (e £ 5(10) £ 5(0 15 + ei£^ (01) £^ (10) ) £^ (0 2) 
+ (e - e 3 - ei 2 )£^ (li0) £| 01) + 3ei 2 £ S([U) £| (li0) £f ((M) 
+3(e 3 - ei 2 )£| (01) £^ (lt0) £^ (0il) + (ei - 2e 3 + ei 2 )£| oi) £ ?(10) j T 
SI 4 r -i 

+ 24 | £ «(0,4) +4£ «(0.1)' C ?(0.3) +3 4(0.2) + 6£ I(0,1) £ «(0.2) + 4(0,1) } T 

+ O s (A,ft). (30) 

In this expression we can see the way in which the parameters {£,4} [real constants 
subject to the conditions 128fl ] describe the arbitrariness we have in the reconstruction 
of the Taylor expansion of a two-parameter family of diffeomorphisms in terms of Lie 
derivative operators. 

To finish this section we will show how to recover the one-parameter case from 
the two-parameter case. The case when one of A or Q vanishes is trivial and it can 
be recovered from the above expressions just setting either A = or f2 = 0. Let then 
consider the only other case, which arises when the two parameters A and f2 are not 
longer independent, e.g. f2 = f2(A). Then the specific way of recovering the single 
parameter (A) case will depend on the specific function f2(A). Here we illustrate the 
simplest case of a linear relation Q — aX (a ^ 0). To arrive to the one-parameter 
expansion (see 1241 125|L we need to study the consequences of the dependence 
between the two parameters. We can do that by looking at the definitions of the 
operators £( p . g ) in equations (|14I27|I . and at their association with Lie derivative 
operators, described by equation (|29|l . The result is a set of relations between the 
vector fields C(p,q) that can be summarized in the following relation: 

= ^£(p+9,o) ■ (31) 

Then, rescaling the vector fields £(p,o) in the following way: £(p.o) — > (p + l) _1 £( Pj o)) 
we arrive at exactly the same expansions as in the one-parameter case [2 fM\ ESI ■ I n 
the case in which both parameters A and fl have a specific physical meaning, it will be 
the physics that will impose the functional dependence fl( A) in particular sub-cases. 



3. Gauges in perturbation theory and the two parameter case 

Let us consider for the moment a spacetime {g D 1 Mo} which we call the background, 
and a physical spacetime {g,A4} which we attempt to describe as a perturbation of 
{g h , A^o}-1f In relativistic perturbation theory we are used to write expressions of the 
form 

9v» ( x ) = ( x ) + {%) , (32) 

relating a perturbed tensor field such as the metric with the background value of the 
same field and with the perturbation. In doing this, we are implicitly assigning a 
correspondence between points of the perturbed and the background spacetimes |26j . 
Indeed through Q32p. which is a relation between the images of the fields in H m 

% As manifolds A4q and M. are the same; for generality we assume that they are m-dimensional. 
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Figure 1. By choosing the coordinates on A4q and M. in such a way that 
(Z o ip)f = x^ , a curve in A^o> and its (£— transformed in Ai have the same 
representation in IR m . Therefore, the components of the tangent vectors V 
and (p*V at the points ip(p) and p are the same: (ip'V)^ (x) = (ip*V)x >1 | = 
V{X o \ v(q) = V z» \ v(q) = V» (x). 

rather than between the fields themselves on the respective manifolds Ai and Aio, 
we are saying that there is a unique point x in JR m that is at the same time the 
image of two points: one (say q) in A4q and one (o) in M. This correspondence is 
what is usually called a gauge choice in the context of perturbation theory. Clearly, 
this is more than the usual assignment of coordinate labels to points of a single 
spacetime [57 . Furthermore, the correspondence established by relations such as (|32|l 
is not per se unique, but rather (|32fl typically defines a set of gauges, unless certain 
specific restrictions are satisfied by the fields involved (e.g. , some metric components 
vanish). Leaving this problem aside, i.e. supposing that the gauge has been somehow 
completely fixed, let us look more precisely at the implications of adopting the 
geometrical description illustrated in Figure ^ 

If we call X the chart on VWo and Z the chart on A4 we see that if we choose 
z(o) = x(q), i.e. the correspondence between points of {g h ,Mo} and {g, M.} implicit 
in 13211 . we are implicitly defining a map tp between A^o and M, such that ip = Z~ 1 oX. 
Thus from the geometrical point of view a gauge choice is an identification of points of 
A^o and Ai. Therefore, we could as well start directly assigning the point identification 
map tp first, calling tp itself a gauge, and defining coordinates adapted to it later. 
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This turns out to be a simpler way of proceeding in order to derive the gauge 
transformations in the following section. 

Let us follow this idea in the specific case of two parameters introducing, in the 
spirit of |23ll21 l2l*ll25| . an (m + 2)-dimensional manifold TV, foliated by m-dimensional 
submanifolds diffeomorphic to TW, so that TV — M. x H 2 . We shall label each copy of 
M. by the corresponding value of the parameters A, ft. The manifold TV has a natural 
differentiable structure which is the direct product of those of A4 and H 2 . We can 
then choose charts on TV in which a; M (// = 0, 1, . . . m — 1) are coordinates of each leaf 
Mam and x m = A, x m+1 = ft. 

In this construction we are assuming that the perturbed spacetimes have the 
same manifold as the background one. With this we are not allowing the possibility of 
addressing questions like how can perturbations affect the differential structure of the 
background spacetime. These issues would require the use of a much more complicated 
mathematical apparatus, in particular the notion of limits of spacetimes introduced 
by Geroch |SS] in order to define the background manifold as a limit A, ft — > of a 
family of perturbed manifolds M\^. Instead, we consider perturbations as fields living 
on the background (as in [23 121 E3I EH]), a standard approach in which these issues 
do not appear. More sofisticated structures of the extended manifold TV have been 
considered in |22j . in an attempt to give the background metric the status of a gauge 
invariant quantity. 

Coming back to our construction, if a tensor T\,q is given on each TVIa.q, we 
have that a tensor field T is automatically defined on TV by the relation T(p, A, ft) := 
T\ t n(p), with p € A4x t n- + In particular, on each AA\si one has a metric g\ t n and a 
set of matter fields r^n, satisfying the set of field equations 

£[g\,n,T\,n} = 0. (33) 

Correspondingly, the fields g and r arc defined on TV. 

We now want to define the perturbation in any tensor T, therefore we must find 
a way to compare T\ t Q with To: this requires a prescription for identifying points of 
■A4x,n with those of A4q. This is easily accomplished by assigning a diffeomorphism 
(f\,n ■ TV — > TV such that p\,n\ Mo ■ TWo — > M\^n- Clearly, <^x,o can be regarded as the 
member of a two-parameter group of diffeomorphisms p on TV, corresponding to the 
values of A, ft of the group parameter. Therefore, we could equally well give the vector 
fields ^i], V C that generate tp. In the chart introduced above, v n m = 1, f^rn+i _ 
v ( m = 0, f (m+i _ ^ but, except for these conditions, v rj : V C remain arbitrary. For 
convenience, we shall also refer to such a pair of vector fields as a gauge. It is always 
possible to take the chart above defined such that v rf = ^Q 1 = 0. So, in this chart, 
point of different submanifolds A4x,n connected by the diffeomorphism ip have the 
same A4 -coordinates x°, . . . ,x m_1 , and differ only by the value of the coordinates 
A, ft. We call such a chart "adapted to the gauge cp" : this is what is always used in 
practice. 

The perturbation in T can now be defined simply as 

^Kn--=<PtnT\ Mo -To. (34) 

+ Tensor fields on Af constructed in this way are "tangent" to A4, in the sense that their components 
m and m + 1 in the charts we have defined vanish identically. 
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The first term on the right-hand side of (|34[1 can be Taylor-expanded using I|1U|) to 
get 



where 



g(k,k') T 



E 



k,k'=0 



Xk ^ k 5 (k,k') T 
k\k'\ ^ 



Qk+k' 



dX k dn k 



£k pk rp 



\=o,n=o,M 



M 



which defines the perturbation of order (fc, k') of T (notice that 5^ T = To) 
worth noticing AqT^q and 8^' k are defined on .Mo; this formalizes the statement 
one commonly finds in the literature that "perturbations are fields living in the 
background". It is important to appreciate that the parameters A, ft labelling the 
various spacetime models also serve to perform the expansion 1)35(1 . and therefore 
determine what one means by "perturbations of order (fc,/V)". 



(35) 



(36) 



It is 



4. Gauge invariance and gauge transformations 

Let us now suppose that two gauges tp and -0, described by pairs of vector fields 
( Vr l: V C) and (^V: ^C) respectively, are defined on TV, such that in the chart discussed 
above* 



( p'n m = 




= l, 




ip £rn 


= o, 






= o, 




ip 


= 1. 



(37) 



Correspondingly, the integral curves of V C) and ^C) define two two- 
parameter groups of diffcomorphisms <p and ip on W, that connect any two leaves 
of the foliation. Thus, (^rj, v () and (^r), are everywhere transverse to M.\,n and 
points lying on the same integral surface of either of the two are to be regarded as the 
same point within the respective gauge: <p and ip are both point identification maps, 
i.e. two different gauge choices. 

The pairs of vector fields ( v n, v £) and (^77, can both be used to pull back a 
generic tensor T and therefore to construct two other tensor fields p* x n T and ip x n T, 
for any given value of (A, f2). In particular, on A4q we now have three tensor fields, 
i.e. To and 

T \,n fXn^Mo ' T tn : = 1>l,nT\ Mo ■ (38) 
Since ip and ip represent gauge choices for mapping a perturbed manifold A4\ t n into 
the unperturbed one Aio, T xn and T xn are the representations, in Mo, of the 
perturbed tensor according to the two gauges. We can write, using (|3^|) - (|3^|l and 
the expansion ifTHjl . 

00 \kctk f 00 \kctk' 

T tn - E w^' fc ' )T = E m £ ^ t = n + A ° T ^ ' (39) 

oc \kr\k f 00 \kok f 

T U = E W^' fc ' )T = E W £ V"'c T = To + AtT x , n (40) 

fe=0 k,k'=0 
* In general, if the chart is adapted to the gauge ip, i.e. v rj M = = 0, it is not adapted to the 
gauge tp, so ^rf ^ 0, + 0. 
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where 8\p' T, 5^ T are the perturbations (|36[) in the gauges p and ip respectively, 
i.e. 



s (k,k') T= C\CK(T 

?(k,k ) r p pk nh 1 rp 

0^ 1 — Ls^L,^! 



(41) 

M 

(42) 

M 



4.1. Gauge invariance 

If T£q — for any pair of gauges ip and ip, we say that T is totally gauge 

invariant. This is a very strong condition, because then (|39l) and 14U|) imply that 
^(k,k ) j, _ ^(fe,fe ) j,^ ^ g au g eg ^ an( j ^ an( j £ or an y ^ j n an y practical case, 

however, one is interested in perturbations up to a fixed order. It is thus convenient 
to weaken the definition above, saying that T is gauge invariant up to order (n, n') iff 
for any two gauges ip and ip 

§{k,h') T = s (k,k') T V (fc,/c') with k < n, k' <n' . (43) 
We have that a tensor field T is gauge invariant to order [n, nf) iff in a given 

c(k k r ) 

gauge cp we have that C^o^ ' T — 0, for any vector field £ defined on .M and for any 
(k,k') < (n,n')< 

To prove this statement, let us first show that it is true for (n,n') — (1,0). In 
fact, if 5p'°^T — S^'^T for two arbitrary gauges tp,ip, we have £^_^T|x = 0. 
But since ip and ip are arbitrary gauges, it follows that v rj — *Vj is an arbitrary 
field f, and f" = £ m+1 = because ^r?™ = = 1, V™ +1 = V' +1 = 0, so 
£ is tangent to In the same way one proves the statement for (n,n') = (0,1). 
Now let us suppose that the statement is true for some (n, n'). Then, if one also 
has 5 { v n+iy) T\ Mo = 8^ +1 ' n,) T\ Mo , it follows that C^^sft'^T = 0, while if 

8tj?' n +1 ^T\mo = ™ +1 '^1a4 ; it follows that C^_^5^" n = 0, and we establish 
the result by induction over (n, n'). 

As a consequence, T is gauge invariant to order (n, n') iff To and all its 
perturbations of order lower than (n, n') are, in any gauge, either vanishing or constant 
scalars, or a combination of Kronecker deltas with constant coefficients. Thus, this 
generalizes to an arbitrary order (n, n') and to the two-parameter case the results 
of [2 E3 1231' Further, it then follows that T is totally gauge invariant iff it is a 
combination of Kronecker deltas with coefficients depending only on A, ft. 



4-2. Gauge transformations 

If a tensor T is not gauge invariant, it is important to know how its representation on 
A^o changes under a gauge transformation. To this purpose it is natural to introduce, 
for each value of (A, fi) e H 2 , the diffeomorphism <J>A,n : Mo — ► Mo defined by 

$A,n := <P-x,si ° i>\,n ■ (44) 

Given that from the geometrical point of view adopted here ip and tp are two gauges, 
$ represents the gauge transformation. The action of $a,o is illustrated in Fig. [3 
We must stress that $ : Mo x IR 2 — > Mo thus defined, is not a two-parameter 
group of diffeomorphisms in Mo- In fact, °®\ 2 ,n 2 ®\ 1 +\ 2 ,n 1 +n 2 , essentially 

because the fields ( v n, v () and (jrj, have, in general, a non-vanishing commutator. 
However, it can be Taylor expanded, using the results of section (|2.2fl . 
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Figure 2. The action of a gauge transformation <E>^ q, represented on the 
background spacetime Mo- The gauge ip is a two-dimensional submanifold 
embedded in M = M X IR 2 , represented by the surface connecting the points 
p, pi, p2, and o. Similarly, the gauge (p is represented by the surface connecting 
the points q, qi, 92, and o. 



The tensor fields and T^q, defined on A4q by the gauges ip and ip, are 

connected by the linear map Q* x n : 

T \,n = ^tn T \ Mo = (iPx,a c P-x,-(i<Px,a T )\ Mo 

= ^nK n r)U=W. (45) 

Thus, the gauge transformation to an arbitrary order (n, n') is given by the Taylor 
expansion of the pull-back <I>^ n T, whose terms are explicitly given in section I^T^l Up 
to fourth order, we have explicitly from (|30l) 

Tt,n = Tin + + 



x,n 

A_ 2 r . ... 1 ... n 2 

~2 



{ A ( 2,o) + £|i. ) } T tn + ~2 { A ( o,2) + £f ( o,i) } T A^f 

Xn + £ oA(i,o)A(o,i) + £ l^(0,i)4 (I ,o) } T A 

A 3 r 1 

~Q (A(3,0) + 3 A(1,0)A(2,0) + ^(1,0) J 



nip 

x,n 



{-^5(2,1) "t" ^^(1,0)^(1,1) ^~ ^€(0,1) ^?(2,0) "I" ^ e 2^f( li0 ) ^£(0,1) ^£(1,0 
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+ (ei - e 2 )C ko l) C\ 1Q) + (e - e 2 )£f (1 0) £f (0jl) } 



rptfi 

1 x,n 



Xfl 2 

9 {^£(1.2) "I" (o,i)^£(i,i) ^£(1.0) ^£(0,2) + ^ e 3^--C(o,i) ^£(1.0) ^£(0,1 



rptfi 

' J A,n 



+ (e - es) A(i,o) A: (0 ,i) + ( £ i " e 3) £ ? ( o,i) ^a.o) } 

. 3 1 

+ "g" \A(0,3) + 3 ^?(0.1)A(0,2) + 4(0,1) J 

+ \ £ ?(4,0) + 4 A (1 ,0) £ S(3,0) + 3 4<2,0) + 6£ |l,0) A=(2,0) + 4(1,0) / ^ 



{^"£(3,1) + 3 -^C(l,0) ^£(2,1) ^ ^£(0,1) ^£(3,0) "I" 3e 4-^-£(2,0) ^£(1,1) 



A 3 » 
"IT 

3e 5^£(l,l)A(2,0) + 3 4(l,0)4l.l) + 3 ( e O^£(l,0)A(0,l) + e l^£(0,l)^£(l,0)) ^£(2,0) 

3 i o, /- z' z-2 

0) 



+ (ei - £2 - e 6 )£^ (Dil) £ fi(10) + 3e 6 £f (li0) £f (0il) £f (1 
+3(e 2 - e 6 )£| (lj0) £ 5(01) £f (1 _ D) + (e - 2e 2 + e 6 )^| (1 0) £f (0il) } T ; 



A. 12 



A 2 ft 2 



| ^£(2, 2) "I" ^£(1,0) "^£(1,2) + ^£(0,1) ^£(2,1) + ^4(1, 



f 

a, a 



+ £ 7^ (2 ,0)-^£(0,2) + e 8-^£(0.2) ^£(2,0) + 4(1,0) ^£(0,2) 4(0, 1) ^(^.O) 

+ 4 ( e o££ ( i l0 )££(o,i) + £ l^£(o,i)' C £(i,o)) ^(i,i) 

-(e 3 + e 2 - ei + e9)4(o,x) 4 ( i,o) - ( £ 3 + £2 - eo - e 9 )4 (1 0) 4 (0 a) 
+ 2(e 3 + e 2 - eoei - e 9)£f( liD) ££ ( o,i)££(i,o)££(o,i) 
+ 2(e 3 + e 2 - eoei + £9)^(0,1) A(i,a)A(o,i)A(i,o) 
-2(e 3 - e ei)£f (liO) £| Ojl) ^ (I(O) - 2(e 2 - e ei)£^ 0jl) £| (i 0) A ( o,i) } T ; 
AO 3 

^ g - {^--£(1.3) + 3 -^£(o ,1)^(1 , 2 ) + ^£(i,o) ^£(0,3) + 3eio^ (0 ,2)^-£(i ,i> 

+ 3e nA(l,l)A(0,l| + 3 4(0,1)A(1,1) + 3 ( £ 04(l,O)^(0,l) + e l^£(0,l)^£(l,0)) ^£(0,2) 

+ (e - e 3 - ei 2 )£f (1)0) £| (01) + 3ei 2 £ ?(M) £ ?(1D) 
+3(e 3 - e 12 )£| 0ii) ^ (10) £ €(01) + (ei - 2e 3 + e 12 )£ 3 (oa) 0) } T A ^ 
fi 4 r 1 

+ ^4 ( £ £(0,4, + 4£ £(0,1) £ £<0,3) + 3 4(0,2> + 6 4(0,l) /: £(0,2) + 4(0,1) } 

+ 5 (\,n), (46) 

where the £(j>, g ) are now the generators of the gauge transformation $A,fi- 

We can now relate the perturbations in the two gauges. To order (n, n') with 
n + n' < 4, these relations can be derived by substituting (|3"5|l . (|4tj|) in i|4T)|) : 

(5 (M) T _ ( ja,o) T = jC ^ o)To) (47) 
(5 y,i) T _ (5 (o,i) T = A(0ii)ro! (48) 
6 (j,o) T _ 6 m T = 2A(i _ o) j(i.o) T + 1 £?(2 o> + } To , (49) 

+ {^-£(1,1) + e o££(i,o)-^£(o,i) + £ i-^£(o,i)-^£(i,o) } To , (50) 
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5f 2) T - S^T = 2C^ r /^T + { A (0 , 2) + 4 0il) } T , 

6 (m T _ 5 im T = 3£?(i o/ (2,o )T + 3 1 £?(2 o) + ^ | 5 (!,o) T 

(3,0) ,0) ^5(2,0) +4, 1,0) 

^(2,D T _ ^(2, 1)T = 2A(i q)(5 (M) T + A(o i) ^,0) T + { ^ Q) + £ 2 i } S (0S )T 

+ 2 {A ( i,D + e o£e (1>0) A(o,d + e iA ( o, i)A(i,o) } <^ 1,0)r 

+ {-^C(2.i) +2^(1,0)^(1,1) + Ai(o,i) A;<2,o) + 2 e 2£^( lj0 )£^( 0il )£{ (li0 ) 

+ (ei - £2)^(0,1)4(1.0) + ( £ ° - e2 ) £ |i,o) £ ?<o,i) } T o . 

^(i,2) T _ d ^ T = 2A(o ^(M) T + A(io) ^) T + { 2) + ^ | ^,o) T 

+ 2 {A(i,d + e oA ( i,o) A ( o,i) + £ i £ c ( o, 1)^(1,0) } <^°' 1)r 

+ {^£(1,2) + 2 Aj(0, l)Aj(l,l) + Aj(l,0) Aj(0,2) + 2e 3A;(0,l) ^£(1,0) ^?(0,1) 

+ (e - e3)A(i,o) £ I(o,i) + ( £l _ e 3) £ ?(0,i) £ ?(i,0) } T , 

§ (o,s )T „ §( o,s )T = 3C , (OA) 5^) T + 3 { £?(o 2) + ^ | ^o,i) T 

+ {^5(0,3) + 3 A(0,1)' C ?(0,2) + £ |o,l) } T , 

6^T - S^T = 4C 5{1A) S^T + 6 {c^ 0) + C\ i o) } 6™T 

+ 4 |A(3, 0) + 3£e (1>0 )£C(2,0) +4(1,0) 

+ {^(4.0) + 4 A(1,0)' C ?(3,0) + 3 4(2,0) + 6 4(1,0)' C C(2,0) + 4(1,0) } T ° , 

5 (3,1) T _ 5 (3,1) T = ^(2, 1)T + ^(3,0) T + 3 { £?(2 Q) + £2^ Q) } ,<1,1> T 

+ 3 {4i,d + e oA ( i,o) A ( o,i) + £ i4o,d4i,o) } 4?' 0)r 
+ (a ( 3,o) + 3 A ( i,o) A ( 2,o) + 4(i,o) } s v' 1)T 

+ 3 {Aj(2,l) + 2 Aj(l, 0)^1(1,1) +^1(0,1)^(2,0) + 2e 2-£C(l, 0)Aj(0, l)Aj(l,0) 

+( C1 - e 2 )A (0il) 4 (1|0) + (e - e 2 )4 (1 , 0) A (0il) } 5™T 

+ {^£(3,1) + 3 Aj(l, 0)Aj(2,l) + Aj(0,l) Aj(3,0) + 3C 4£C(2,0) ^4(1,1) + 3e 5^£(l, 1)^4(2,0) 
+ 3 4(1,0) + 3 ( £ 0-CC(l, 0)^(0,1) + e lA(»,i)A(i,o))A(2,o) 

+ (ei - £ 2 - ee)^,!) 4(1,0) + 3e 6A(i,o)A ( o,i) 4(1,0) 

+3(e 2 - e 6 )4 (li0) A(o, uAa, o) + ( £ o - 2e 2 + e 6 )4 (lj0) A(o,i) } T o > 
^(2,2) T _ ^ 2 , 2)T = 2£?(i ^(i,2) T + 2£ C(M) a< a . 1 >T 

+ {A(2,o, + 4(1,0) } ^ 2) T + {^(0,2, + 4<o,i) } 
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+ 4 + e 0£?(l,o)£? (O ,l) + £ lA(O,l)A(l,0)} V' T 

+ 2 {£5(2,1) +^5(1,0)^(1,1) +^5(0,1)^-5(2,0) + ^ e 2^5(l, 0)Aj(0, 1)^5(1,0) 

- 6 9 )A (0il) 4 M) + ( £0 - 6 2 )£|, 0) £ 5(0 , 1) } S^T 

+ 2 \£5(i,2) + ^^5(o,i) ^5(1,1) + -^5(1.0) -^5(0,2) + 2 £ 3^ (0 ,1)^5(1 ,0)^5(0 ,i) 
+ ( £o - e 3 )A (1 , 0) £| (0>1) + ( Cl - e 3 )4 (0>1) A (1 , , } ^' 0) T 

+ {^5(2.2) + 2 A(1,0)A(1,2) + 2 ^5(0,1)A(2,1) + 2 4(1,1) 

+ € ~!^'i{2,0) •^'5(0,2) + £ 8£C(o, 2)^5(2,0) + 4(1,0) ££(0,2) + 4( , 1) Aj(2 ,0) 

+ 4 (eoA ( i, 0)^5(0,1) + e i^5<o, 1)^5(1,0)) 

~(C3 + £2 - Cl + £9)4(0,1)4(1,0) ~ ( £ 3 + «2 - £0 - £9)4(1,0)4,0,1) 

+ 2(e 3 + e 2 - £o£i - e 9)£? ( i,o)£5(o,i)£5(i, o)££(o,i) 
+ 2(e 3 + e 2 - e £i + e 9)£? ( o,i)£5(i,o)£5(o, i)A(i,o) 

-2(e 3 - e ei)£5(i.o)4(o,i) /: 5(i,o) _ 2 ( e 2 ~ e o e i)£5(o,i)4(i,o ) ' C 5(o,i) } T o i ( 58 ) 
5 (i,3) T _ §( i, 3)T = 3A(o ^(x, 2)T + A(i d/ (o,3) t + 3 | A(q 2) + ^ I 5 (i,i) T 

+ 3 {£5(i,d + e o£5(i.o ) ' C 5(o,i) + £ i £ 5 ( o,i) /: 5(i,o)} ^°' 2)r 
+ { £ 5(o.3) + 3A(o,i) A(o,2) + 4(o,i) } S v'° )t 

+ 3 {-^5(i,2) + 2 -^5(o, 1)^5(1,1) + ^-5(1,0) ^-5(0,2) + 2e 3£c (0 ,i)£C(i, 0)^5(0,1) 

+ ( £0 - 63)^(1.0,4(0,1) + ( £ 1 ~ e 3)4 ( o,l,A(i,o, } 

+ {-^5(1.3) + 3-^5(0, 1)^5(1, 2) + ^5(1, 0)^5(0, 3) + 3eio£5 (0 ,2) £5(1,1) + 3e n££(i ( i)££(o ( 2) 

+ 3 4(o,i)£5(i,i) + 3 ( e O-C£ ( i,o)££(o,i) + e l£5 ( o,i) £5(i,o)) ^5(0, 2) 
+ (e - £3 - ei 2 )A(i,o) 4(o,i) + 3e 12£5<o.i, £5(i.o, 4(o,i) 

+3(e 3 - e 12 )C\ o i) C kl 0) £ C(0 1} + ( Cl - 2e 3 + ei 2 )4 (ol) £ €(M) } T , (59) 
sf^T - S^T = iC^S^T + 6 {A (0 , 2) + 4 (M) } S^T 

+ 4 {£5(o,3) + 3 £5(0,1) £5(0,2) +4(0,1) 

+ {£5(0., + 4A«o,i, £c«o,3, + 3 4,o,2, +64(o,i ) £5 ( o,2, +4,0,1,} To- (60) 

This result is, of course, consistent with the characterization of gauge invariance given 
in subsection 14. II Equations l|4TJl and imply that Ta.q is gauge invariant to the 
order (1,0) or (0, 1) iff C^Tq — 0, for any vector field on .Mo- Equation (|49|l implies 
that T\ t Q is gauge invariant to the order (2, 0) iff To = and C^S^^T = 0, for any 
vector field on Mo, and so on for all the orders. 

It is also possible to find the explicit expressions for the generators £,(p.q) °f the 
gauge transformation $ in terms of the gauge vector fields ( v r), v () and (^77, ^£). We 
write here their expressions up to second order: 

£(1,0)=^-^, (61) 



Two-parameter non-linear spacetime perturbations 



18 



m 



x(M ) 



X.' 



Mo 



Z 



<l 



y(q) = x(p) 



x{q) 




A 




P 







Figure 3. The picture shows how two different gauge choices, ip and ip, map two 
different points in the background manifold Mo, namely q and p respectively, to 
the same point o in the physical spacetime M\n- The coordinate version of this 
fact is shown on the left part of picture. We can look at this either from the active 
point of view (the two points q and p are given coordinates in the same chart by 
the map X, namely x^(q) and or from the passive point of view (the point 

q is given coordinates in two different charts by maps X and Y, corresponding to 
two different points of ]R m , x(q) and y(q) = x(p)). Finally, Z assigns coordinates 
to the point o of the physical spacetime M\q. 



4-3. Coordinate transformations 

Up to now, we have built a two-parameter formalism using a geometrical, coordinate- 
free language. However, in order to carry out explicit calculations in a practical case, 
one has to introduce systems of local coordinates. In this respect, all our expressions 
are immediately translated into components simply by using the expression of the 
components of the Lie derivative of a tensor. Nonetheless, much of the literature on 
the subject is written using coordinate systems, and gauge transformations are most 
often represented by the corresponding coordinate transformations. For this reason, 
we devote this subsection to describe how to establish the translation between the 
two languages, giving in particular the explicit transformation of coordinates (further 
details are in 3 ). 



£(0,1) 

£(2,0) 

£(i,i) 

£(0,2) 



(62) 
(63) 
(64) 
(65) 



[Vrj], 
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Let us consider the situation described in Fig. |3J We have considered two gauge 
choices, represented by the groups of diffeomorphisms <p and tp, under which the point o 
in the physical manifold M.\ ^ corresponds to two different points in the background 
manifold Mo, namely q = ^n(°) anc ^ P ~ ^\l~i( )- The transformation relating 
these two gauge choices is described by the two-parameter family of diffeomorphisms 
$A,n = (f\ q V'A.n, so that <I>A,n(p) = q. This gauge transformation maps a tensor 
field T on q e A4 to the tensor' field ($*T)(p) = $*(T(g)) on p G A4 - 

Now, let us consider a chart (U,X) on an open subset W of Mo- The gauges y>A,fi 
and ipx.fi define two maps from M.\,n to IR m : 

Xo^: M x , n — H m ^o^j,: A<1a,q — 

o | — > z(<?( ))> o | — ► a;(p(o)). 



Then, we can look at the gauge transformation $ in two different ways: from the 
active point of view or from the passive point of view. In the first case, one considers a 
diffeomorphism which changes the point on the background Mo- To these points one 
associates different values of the coordinates in the chart (U,X). So the coordinate 
change is given by 

x»(p) — x»(q) (67) 

or, defining the pull-back of x as x M (p) := x M (<I>(p)), 

a^(p) — ► x"(p). (68) 

If we instead consider the passive point of view, we need to introduce a new chart 
(U',y): 

Y:=Xo , (69) 

such that the two sets of coordinates are related by 

y»(q) = xHp) , (70) 

so we can say that the gauge transformation does not change the point on Mo, but 
it changes the chart from (U,X) to iU' ,Y), i.e. the labels of the point of Mo- The 
coordinate transformation is then 

x*(q) — y»(q). (71) 

Now, let us consider the transformation of a vector field V. From the active point of 
view, the components of V in the chart ilA,X), V^, are related with the ones of the 
transformed vector field V, V 1 *, by 

V» = (X*Vr = (X*$l Q Vr. (72) 

From the passive point of view, we can use the properties relating the pull-back and 
push-forward maps associated with diffeomorphisms: 

X,** x>n V = X^l n V = Y*V , (73) 

so we get the well known result that the components of the transformed vector V in 
the coordinate system X are defined in terms of the components of the vector V in 
the new coordinate system Y: 



r". ,■(/.). = (Y- i r. (/ ))" = r''-( // (^.= ( g£ 



V v {x{q)). (74) 

x(q) 
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In order to write down explicit expressions, we will apply the expansion of the 
pull-back of $* [See equation i|30|) ] to the coordinate functions . Then, the active 
coordinate transformation is given by 

x K(p) =x"(q) = (#V)(p) 

A 2 / \ ft 2 / \ 

+ ~2 V^( 2 >°) "'"^( 1 . )^(i>o),fJ + ~2 v(°' 2 ) ~*~ ^(0.1)^(0,1),^ 

+ "g~ (^3,0) + 3 £(l,0)£f2,0), I / + £fl,0)ffl,0) 1 />£(l,0),i/) 

A 2 ft / 

"I 2~ V( 2 '!) + 2 £(l:0)£(l,l),y + S(0, l)£(2,0),i/ + 2e2 £(l,0)£(0,l),p£(l,0),^ 

+ (ei ~ e 2)C('o,i)C(i,o),p^(i,o)^ + ( £ o ~ e 2)£(i,o)ffi,o),pf(o,i),i/) 

(^(1,2) + 2 £(0>l)£fl,lV + £(1>0)£(0,2),K + 2e3 ^(0,l)^C l ..O),^(0,l),i/ 



Aft 2 

+ — 

+ (eo - e3)£f 1)0) £(V),p£(o,i)," + ( £l ~ £3 )£(V)£(V),p£(i>oV 

Q3 , ^ > 

+ ~g~ (^(0,3) + 3 £(0, l)£(0,2),i/ +£(0,l)£(0. 1 ).p£(0,l),i'J 

A 4 / 

+ ^^(i,o)^ri>o)>/>^(2,o),i/ +^fi,o)^(i,o),o-^r i .o).p^ri,o),i/)) 

A 3 ft / 

^ ^~ (£(3,1) + 3 ^a.0)^(2,l),i/ + £(0,l)f(3,0),i/ + 3e4 ^r2.0)^(l,l), 1 / 

+ 3e 5^fl, 1)^(2,0),!/ + 3 ffl,0)£Fl,0),p£(l,l),i/ + 3 ( £ °£(1,0)£(0,1),P + e l£(Q,l)£Fl,0),p)£(2,0),K 

+ (ei - £2 - e e)C(Q,i)^(i 1 o) )CT C(i 1 o),pC(i ) o), I / + 3e 6^fi,o)C(o,i), CT ^a,o),p^(i,o),i/ 

+ 3(e 2 - e 6)£(l,0)£(l,0),at(0tf,p€(lfi),v + ( £ - 2e 2 + e 6)^fl,0)^fl > 0) lff ^( , l,0) 1 pCfo ) l), I /) 

A 2 ft 2 / 

^ — {£(2,2) + 2 ^(1.0)^(1,2),(/ + 2 £(0,l)€(2,l),u + 2 £(i,l)£(l,l),i/ 

+ e 7£(2,0)£(0,2),i/ + e8 £fo, 2)£(2,0),k + £(1, 0)£(1> 0),p£(0,2),i/ + £(0,l)£(0,l),p£(2,0),i/ 

+ 4 ( e 0^ lj0 )?(0,l),P + e l£(0,l)£(l,0),p)£(l 

- (£3 + £2 - Cl + £9)^0,1)^0,1)^^1,0)^^1,0),!, 

- (£3 + £2 - £0 - £9)^1,0)^(1,0)^^0,1)^^0,1),!/ 

+ 2(e 3 + £ 2 - £ £l - £9)£(%)£f ,l), CT £(%),p£(0,l)," 

+ 2(e 3 + e 2 - e £i + £9)^0,1)^1,0)^^0,1)^^1,0),!/ 

- 2(e 3 - £o£i)^(i,o)C(o,i) !0 .C(o,i),p^(i,o), y ~ 2 ( £ 2 - £o£i)C(o,i)Cfi,o), (T ^ri,o),/ 9 ^(o,i), ! /) 
Aft 3 / 

H 6~ V^(l= 3 ) + 3 £(0>1)£(1,2),«/ + £fl,0)£(0,3),i/ + 3e 10C(0,2)^(l,l), ! / 

+ 3e llC(i,l)f(o,2),i/ + 3 £(0,l)£(0,l)>p£(l,l)," + 3 ( e °£(l,0)£(0,l),P + £ l£(0, l)£(l,0),p)£(0,2), !/ 
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+ (eo e 3 e i2)^(i,o)^(o,i),<T^ro>i),p^ro,i)^ + 3ei2 ^ro,i)^(i,o),-T^( I o,i),p^(o,i)^ 




where the vector fields & p q \ and their derivatives are evaluated in x(p). This 
expression gives the relation between the coordinates, in the chart (U,X), of the 
two points p and q of Mq. 

On the other hand, the passive coordinate transformation is found by 
inverting (|731) : 

y»(q) := x"(p) = x»(q) - ^ 1>0) (x(p)) - il^ 01) (x(p)) + 2 (X,il) , (76) 

and then by expanding x(p) around x(q). We obtain in this way an expression of the 
form 

y»{q) = x»{q) - X^ m (x(q)) ~ ^ 0A) (x(q)) + 2 (A,tt) , (77) 

which gives the relation between the coordinates of any arbitrary point q s Mo in the 
two charts (U,X) and (W,Y). Such a relation is needed to find the transformation 
of the components of a tensor field, by using , as it is usually done in textbooks 
for first order gauge transformations |291 13(J| . However, in order to determine these 
transformation rules it is much simpler to apply directly the expressions Q47I60JI . 
computing explicitly the Lie derivatives of the tensor field. 

5. Conclusions 

Many astrophysical systems (in particular, oscillating rclativistic rotating stars) can be 
well described by perturbation theory depending on two parameters. A well-founded 
description of two-parameter perturbations can be very useful for such applications, 
specially in order to handle properly perturbations at second order and beyond. For 
example, one may wish to compare results derived in different gauges. 

In this paper we have studied the problem of gauge dependence of non-linear 
perturbations depending on two parameters, considering perturbations of arbitrary 
order in a geometrical perspective, and generalizing the results of the one-parameter 
case j3 E3] to the case of two parameters. We have constructed a geometrical 
framework in which a gauge choice is a two-parameter group of diffeomorphisms, 
while a gauge transformation is a two-parameter family of diffeomorphisms. We have 
shown that any two-parameter family of diffeomorphisms can be expanded in terms 
of Lie derivatives with respect to vectors £^ >. In terms of this expansion, which can 
be deduced order by order, we have derived general expressions for transformations 
of coordinates and tensor perturbations, and the conditions for gauge invariance of 
tensor perturbations. We have computed these expressions up to fourth order in the 
perturbative expansion, i.e. up to terms X k fl k with k + k! = 4. 

The way in which the expansion of a two-parameter family of diffeomorphisms 
was derived in this paper is order by order, constructing derivative operators that can 
be rewritten as Lie derivatives with respect some vector fields. The development of an 
underlying geometrical structure, analogous to the knight diffeomorphisms introduced 
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in the one-parameter case 2 , would be interesting for two reasons: first, in order to 
have a deeper mathematical understanding of the theory, and second, in order to derive 
a close formula, valid at all orders, for gauge transformations and gauge invariance 
conditions. The present paper has been devoted to the derivation of the useful formulae 
for practical applications. In particular, our expressions will be useful to compare 
results derived in different gauges, and can form the basis for the construction of a 
gauge invariant theory of two-parameter systems in the line of works done for the 
one-parameter case like as for example ED E3 EH EH EH ■ We leave the development 
of a more formal framework for future work. 



Appendix A. Proof of the statement (|29B 

The aim of this Appendix is to give a proof of a theorem that allows us to make the 
statement contained in equation i|29|) . 

Theorem: "Let £ be a derivative operator acting on the set of all the tensor fields 
defined on a differentiable manifold Ai and satisfying the following conditions: (i) 
It is linear and satisfies the Leibniz rule; (ii) it is tensor-type preserving; (iii) it 
commutes with every contraction of a tensor field; and (iv) it commutes with the 
exterior differentiation d. Then, there exists a vector field £ such that C is equivalent 
to the Lie derivative operator with respect to £, that is, £g." 

First of all, notice that the operators introduced in equations (|14I27|) satisfy the 
conditions of the theorem. In particular, properties (iii) and (iv) follow from the fact 
that commutes with contractions and the exterior derivative (see ES]). For more 
details on this question see, e.g., pT7] . 

The proof of the theorem is as follows: When acting on functions, C defines a 
vector field £ through the relation 

£/=:£(/), V/e^M), (A.l) 

where T(M) denotes the algebra of C°° functions on M.. What we want to prove is 
that on an arbitrary tensor field T, 

CT=£^T. (A.2) 

Clearly (|A.2(I holds for an arbitrary tensor field T iff it holds for an arbitrary vector 
field V. For the latter equation (| A. 2|> is equivalent to the following expression 

CV = %V\. (A.3) 

Applying this to any function / we obtain 

(£V)(f)=Z[V(f)]-V[af)}, V / G F{M) . (A.4) 

Therefore, to prove (|A.2|I is equivalent to prove (|A.4|I . To this end, let us consider the 
action of the operator £ on the function V(f). Using IjA.lfl we have 

C[V(f)]=aV(f)}. (A.5) 

On the other hand, using the properties (i)-(iv) of C we have 

C[V(f)] = C(df(V)) = £{C(df ® V)] = C[d(£f) ®V + df®CV] 

= d(Cf)V + df(CV) = V(Cf) + (CV)(f) . (A.6) 

Then, this in combination with (|A. 5|) . and using (|A.1|) . leads to equation (|A.4(1 . which 
is what we wanted to prove. 
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